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We consider the effect of mobile ions on the applied potential needed to reorient a lamellar system 
of two different materials placed between two planar electrodes. The reorientation occurs from a 
configuration parallel to the electrodes favored by surface interactions to an orientation perpendicular 
to the electrodes favored by the electric field. The system consists of alternating A and B layers 
with different dielectric constants. The mobile ions are assumed to be insoluble in the B layers and 
hence confined to the A layers. We find that the ions reduce the needed voltage most strongly when 
they are constrained such that each A lamella is electrically neutral. In this case, a macroscopic 
separation of charge and its concomitant lowering of free energy, is attained only in the perpendicular 
orientation. When the ions are free to move between different A layers, such that charge neutrality 
is only required globally, their effect is smaller and depends upon the preferred surface interaction 
of the two materials. Under some conditions, the addition of ions can actually stabilize the parallel 
configuration. Our predictions are relevant to recent experiments conducted on lamellar phases of 
diblock copolymer films with ionic selective impurities. 
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I. INTRODUCTION 

Block copolymers are polymeric materials with specific chain architecture in which several blocks of different chem- 
ical nature are linked covalently. The simplest design is that of a di-block in which two polymeric blocks, A and B, are 
connected. Even in this design the system manifests a rich phase structure of self-assembled periodic arrays of spheres, 
cylinders, or planes, which has many applications such as templates for nanoscopic device*^. One difficulty which 
must be overcome before block copolymer arrays can be exploited is the fact that, when produced as a thin film, the 
array tends to align with the substrate due to preferential interactions between it and one of the polymer blocks. In 
many applications, however, the desired orientation is one in which the ordered structure is aligned perpendicular to 
a substrate. For example, if a cylindrical array is aligned perpendicularly, then the cylinder cores can be etched away 
and replaced by a metal. Further removal of the material surrounding the cores produces an ar ray of wires useful 
for sensors^. One method to bring about this orientation is to subject the film to an electric fielcP^tL^, a technique 
which takes advantage of the difference in dielectric constants between the two polymer blocks. Because the electric 
field needed to produce the desired orientation can be rather high, on the order of 10V//zm, it can happen that the 
material undergoes dielectric breakdown before alignment is achieved. As a way of lowering the electric field needed 
for alignment, Tsori et alP^ suggested that one could introduce free ions into the system. Presumably if the ions were 
preferentially soluble in one of the blocks, the induced dipole moment of that block would become much larger in the 
perpendicular orientation. The increase in polarization would bring about the desired decrease in the electric field 
needed for alignment. 

T his p rocedure has been utilized experimentally and the desired increase in orientation obtained for a given applied 
fiel d 14 * 15 *. In spite of this success, it is not completely clear how the addition of ions brings about the desired result. It 
has been argued that instead of increasing the induced dipole moment, large amounts of added ions introduced in one 
of the two blocks simply change the respective dielectric constant^. A different argument!^ is that the addition of 
large amounts of ions changes the strength of the interaction between the different monomers even in the absence of 
an external field. In support of this latter suggestion, introduction of ions has b een shown to affect the phase behavior 
of block copolymers, behavior which scales with this interaction strengtri- 17 * 18 *. 

To elucidate the origin of free-ion effects in block-copolymer systems, we have considered the effect of free ions 
on a system which is simpler, but which shares its mesoscopic periodic morphology. Our model system is a rigid 
stack of alternating A/B lamellae. The entire stack is confined between two external flat electrodes. The A and B 
lamellae are characterized by two dielectric constants, ka and Kb, respectively. Free ions are then introduced into the 
A layers. It is assumed, arbitrarily, that the positive charges are immobile and distributed uniformly throughout the 
A layers, creating a uniform charged background, while the negative counterions are mobile throughout the A regions. 
An interchange of the role of positive and negative charges will have no effect on the results of our calculation. The 
fraction of ions dissolved in the two kinds of layers is, in general, different, and we consider the case of greatest contrast 
in which there are no charges in the B lamellae. We calculate, within mean-field theory, the system free-energy in the 
presence of a voltage Vq applied across the plates. The free energy depends, of course, upon the relative orientation 
of lamellae and electric field. 

Our results for the critical reorienting electric field depend strongly upon the location of the counterions. If they are 
confined to the A lamellae such that each is electrically neutral, then separation of charge occurs only over distances 
of the order of the thickness of the A lamellae when the lamellae are parallel to the substrate. In contrast, the charges 
can be separated over lengths of the order of the entire film thickness when lamellae are oriented perpendicularly. 
Hence, the latter configuration is favored, and the desired orientation can be brought about with an electric field 
significantly weaker than that needed without ions. However, if the counterions are free to move between different A 
lamellae, then the separation of charge occurs over distances of the order of the sample thickness in both orientations, 
so that the system behavior is not easily predicted. In this case, our results show that there is a smaller reduction of 
the field needed to reorient the sample, and that the magnitude of the reduction depends upon whether the surface 
interactions prefer the A component, whose lamellae contain ions, or the B component. The effect is larger if the 
surface interactions prefer the B component. If they prefer the A component, we find that the introduction of ions 
can even stabilize the parallel orientation. 

II. THE MODEL 

We consider a system, of volume f2, consisting of alternating layers of materials denoted A and B with dielectric 
constants ka and k g , respectively. All layers but the two closest to the two surfaces (chosen arbitrarily to be B layers) 
are of equal thickness A/2 so that the structural periodicity is A. The two layers adjacent to the plate electrodes are 
each of thickness A/4, as they would be in a system of unstrained lamellar block copolymer. The material fills the 
space between two plates, which are parallel to the (x,y) plane, one at z — at which V(x,y,0) = 0, the other at 
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z — d at which V(x, y, d) = Vq. Each plate is of area <S p i a tc- Just as in the experiments, there is a thin insulating layer 
between the plates and the material, keeping the total amount of ions in between the two plates fixed. We assume 
the system contains N + immobile, monovalent, cations distributed with uniform density p + in the A lamellae, and 
zero density within the B lamellae. The counterions are also monovalent and assumed to be mobile with a non-zero 
number density, p_(r), only within the A lamellae; there are no counterions in the B lamellae, and the entire lamellar 
system is electrically neutral. 

The Maxwell equations which govern this electrostatic system with linear dielectric materials are, within the SI 
system, 

V x E(r) = and (1) 
V • (K(r)eoE(r)) = e(p+(r) - p_(r)), (2) 

where e > is the unit of charge. The first equation is satisfied identically by introducing the electrostatic potential 
V(r) such that E(r) = — W(r). With this the remaining Maxwell equation takes the form 



V • («(r)eoW(r)) = -e(p+(r) - M'))- (3) 

In our model system, the number density of positive charges within the A lamellae is constant, p+(r) = p + , and the 
dielectric constant «;(r) = ka, a constant, so that Eq. ([3| becomes the Poisson equation 

V 2 V(r) = —(p + — p-(r)) r in lamellae A. (4) 

Within the B lamellae, the number density of all charges, cations and anions, is taken to vanish, and ft(r) = Kg a 
constant, so that within these regions Eq. ^ becomes the Laplace equation 

V 2 V r (r) = r in lamellae B. (5) 

At an A/B interface at which the dielectric properties of the material changes abruptly, the potential V(r) is continu- 
ous, the parallel component of E is continuous as follows from Eq. 0, and the normal component of the displacement 
field, D = keqE, is continuous as follows from Eq. |2]). The potential is now completely specified as a functional of 
the unknown charge density, 

V(r) = V[p-(r)]. (6) 

Within mean-field theory, the relation between the ensemble-averaged charge density and potential satisfies the 
Boltzmann distribution: 

, , exp[/3eF(r)] . 

p_(r)=p + n VVH \; J rmA, 7 

J A dr expipeV (r)J 

with p = 1/fegT, ks is the Boltzmann constant, T is the absolute temperature, and the integral is taken over all the 
A lamellae. 

Equations Q, (JjjJ, and |7| constitute the self-consistent, Poisson-Boltzmann, equations which must be solved subject 
to the boundary conditions given above. Once the solutions for the charge density and potential are determined, the 
free energy of the system can be obtained. It is convenient to display explicitly the contribution to it from the surface 
fields and to write the total free energy as 

F tot (T, n, V , N+,S A , S b ) - F el + 1A S A + J B S B , (8) 

where Sa is the total area of the two plates in contact with the A region and similarly for Sb', 7a and jb are the 
interfacial free energies per unit area between the plates and A and B regions, respectively; Sa + Sb = 2iS p i a t e - 

For convenience, we remind the reader in Appendix I of the result for the electrostatic energy of the system when the 
plates are held at a fixed potential difference. With the energy in hand, the mean-field expression for the electrostatic 
part of the free energy can be derived in various ways^^H with the result 

^ el (T,n,^ 0) iV + ,<S pla te) = -\e J n{v){VV{v)) 2 dv + e j V{v){p+ - P -{v))dv 

+ k B T [ p_(r)ln(p_(r)/p+)dr. (9) 
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The free energy is proportional to the area of each plate, so it is convenient to consider the free energy per unit 
area. We choose to measure the free energy per area in units of 26qka/ ' /3 2 \e 2 . Hence we introduce dimensionless free 
energies per unit area, 

/3 2 Xe 2 F t t 

ftot{T,p + ,Vo,S A /Spl a te) = ~ c - , (10) 

ZeoKA Opiate 

fa(T, P+ ,V Q ) ee f^^, (11) 

ZCqKa Opiate 

with F e i given by Eq. (JoJ) , so that within mean-field theory, Eq. ^ becomes 

/tot = /el + 7, c 1A + c 7S • (12) 

\ZeoKAj \ Opiate O p i a tc / 

To determine the orientation of the lamellae in the presence of the external field, we determine the free energy for 
the two orientations: lamellae that are perpendicular to the substrate and those that are parallel to it. In the latter 
case, we approximate the situation expected if the system were an unstrained A/B diblock copolymer film. As the 
substrate prefers one of the two blocks, the lamellae next to the two plates consist of the preferred block (assumed 
here to be the B block) and will have a thickness of A/4, while all other A and B lamellae are of thickness A/2. 

We consider two different constraints on the location of the mobile counterions. In the first, we assume that the 
countcrions are not free to move between A lamellae. Their distribution is such that each A lamella is electrically 
neutral. We refer to this as local neutrality. The other possibility we consider is that the countcrions, while found 
only in the A lamellae as before, are distributed among them subject only to the weaker constraint that the system is 
overall neutral. We refer to this as global neutrality. The free energy of the system in which the lamellae are oriented 
parallel to the plates is affected significantly by the difference between these two constraints, while the system of 
perpendicular lamellae is not affected at all. 

Before discussing the solution of the equations, we note that there are a few dimensionless ratios which parametrize 
them. The simplest of these is a dimensionless applied potential, v 

v ee e/3V . (13) 

As typical potentials across the film are a few volts to a few dozens volts, this parameter is large, on the order of 10 2 
for experiments at room temperature. A dimensionless measure, r, of the ionic number density p + is readily defined 



f3e 2 P+ d 2 ^ 1/2 



ka£q 

and is easily seen to be the ratio of the system thickness d to the Debye length, A^ 

d 



(14) 



(15) 



where 

1/2 

(16) 



A / 



3e 2 p + 



Thus, small values of r imply weak screening of the electric field by the counterions, while large values of r imply 
strong screening. The temperature-independent ratio 

r 2 = ep+d 2 ^ p+ 
v K A e V V ' 

depends only on the ratio p+/Vo, which indicates that the behavior in the limit of large externally applied potential 
is the same as for small ionic charge density. 

It is convenient to write the Poisson-Boltzmann equation in dimensionless form. To this end we measure all 
distances in units of d, the distance between the two bounding electrodes, f = r/d, and define the dimensionless 
potential W(x,y,z) = f3eV(x,y, z) so that W(x,y, 0) = 0, W(x, y,l) = v. Then the Poisson-Boltzmann equation 
takes the form 

V 2 W(x, y, z) = -r 2 ( 1 - 9e w{i) \ in lamellae A, (18) 



9 = P - {W =°\ (19) 
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and, of course, Laplace's equation 

V 2 W(x, y,z)=0 in lamellae B. (20) 

We obtain estimates of the parameter r, and therefore the density of free ions from experiments on polystyrene- 
polymethylmethacrylate (PS-PMMA) diblock copolymers where lithium salts are infused in the MMA blocks. It is 
known that Li + ions are associated with the carbonyl groups of the MMA (dielectric constant ka = 6.3). The number 
density of the MMA monomers is about 7 x 10 27 m~ 3 . What is known less well is the number density of Li + ions 
solubilized in the MMA blocks. Wang et alP^ estimate that the largest fraction of Li + is about 0.27, which would 
correspond to a number density p + ~ 2 x 10 27 m~ 3 . If all this charge were mobile, it would yield a nominal Debye 
length Xd = 0.067nm at T — 298 K. A film thickness of 300nm as in Xu et alP^l, would imply a large value of r ps 4500. 
However, such a sub-atomic value for A^> would indicate that continuum theories, such as ours, would be inapplicable. 
In other experiments, the fraction of ions is clearly much smaller. For example, Kohn et al! 2 ^ quote a Debye length, 
\d = 2.3 nm which, for a film of the same thickness, would yield r w 130. Further Tsori et alP^ estimate the the 
fraction of residual Li + ions to MMA to be much smaller than in the work of Wang et alP^, about 3xl0 -5 . This 
means that p + m 2 x 10 23 m~ 3 , corresponding to a Debye length of 6.7 nm. With the same film thickness of 300nm, 
r ~ 45. For films only a few lamellae in thickness, r would be even smaller. Thus the parameter r is expected to vary 
between values greater than, but of order of unity to those as large as the order of 10 3 . For the smaller values in this 
interval, one expects the effect of the free ions to be but a perturbation on the ion-free results, while this is certainly 
not expected to be the case for the larger values. Thus, experiments span a range over which the effect of the ions is 
expected to vary from negligible to important. 



III. RESULTS 



A. The Parallel Orientation 



We first consider the case in which the lamellae are parallel to the plates. The total amounts of A and B material 
are equal, and the B layers have no charges. Were there also no charges in the A layers, the electric fields would be 
E A = 2Vo/[(l + K)d\ and E B = kEa with k = ka/ k b- The potential would increase linearly across the lamellae. The 
magnitudes of the dimensionless electric fields, —W'{z) would be W' A = 2v/(l + k), W' b = kW' a . We consider here 
v = 10, and k = 2. (For PMMA, k a = 6.3 and for PS, k b ^ 2.52, so that « = 2.5.) Therefore, the magnitude of 
the dimensionless electric fields would be W' A = 20/3, W' B = 40/3. Given a system in which there are free ions, we 
must solve the Poisson-Boltzmann and Laplace equations, Eqs. (jl8h and (20 1. This is done numerically by means of a 



procedure described in Appendix II. The solution yields the potential, electric field, and charge density as a function 
of position z. 



1, Local Neutrality 



Consider the case in which each A lamella is electrically neutral. It is clear that in this case a dipole moment scaling 
with the distance between plates can only be obtained when the lamellae are oriented perpendicular to the plates. 
Thus, the applied field will favor this orientation, and we expect the critical field needed to bring it about will be 
reduced by the presence of the free ions. 

We examine a system in which the wavelength of the lamellar structure is A = d/4 because the effect of the ions is 
more easily seen in thin films. There are four layers of A, each of thickness d/8, three layers of B, each of thickness 
d/8 and two layers of B next to the two plates, each of thickness d/16. Were there no ions, the dimensionless potential 
W would be as shown in Fig. [T] The presence of free ions does not noticeably affect the potential, even for values 
of the charge density corresponding to r 2 = 100. Presumably this is because the charge separation can never exceed 
A/2 and so is never macroscopic. The effect on the electric field ~ W , is small, but is discernible for a large charge 
density r 2 = 100, and is shown in Fig. [2} The non-zero charge density, ~ W" , can also be discerned, and is shown in 
Fig. [3j One sees that the variation of the charge density across each A lamella is essentially linear for small charge 
densities, being positive on one side of the lamellae and negative on the other. This is a result of the imposed charge 
neutrality of each lamella. For large charge densities, deviation from this linear behavior is expected. 
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2. Global Neutrality 

Next consider the case in which the A layers are not constrained to be locally neutral; the system is only subject 
to the weaker constraint that the system be globally neutral. While the perpendicular orientation of lamellae is 
unaffected by the difference between these constraints, the parallel orientation is greatly affected. Each A lamella 
will now exhibit a non-zero charge density profile that is a function of the distance to the plates. Thus, there will be 
a dipole moment that scales with system size in the parallel orientation as well as in the perpendicular orientation. 
As a consequence, the latter has no obvious advantage over the former, and it is unclear whether there will be any 
reduction in the field needed to bring about a reorientation of lamellae. 

Results for the potential, electric field and charge density in the parallel orientation are shown in Figs.[4j[5] and[6j 
respectively, for various charge densities. For the small density of charges, r 2 — 0.01, the potential is imperceptibly 
perturbed from the case of no ions, shown in Fig. [lj and the electric field hardly differs from the values 20/3 and 40/3 
which would be obtained were there no ions. We also note that the electric field when r 2 = 100 is not symmetric 
about the midpoint of the sample, z/d = 0.5, because there is no symmetry on interchange of positive ions, which are 
fixed, and negative counterions, which are free to move. 

As the density of free ions is increased, the charges are free to migrate towards the bounding electrodes and to 
concentrate at the surfaces of the lamellae. The charge densities are much larger than in the locally neutral case. 
(Compare the scales of the charge density in Fig. [6] with those in the locally neutral case, Fig. [3]). One notes from 
Fig. [6 that, for large charge densities, r 2 , the mobile negative charges are essentially depleted from the lamellae close 
to the negatively charged plate leaving behind a positive charge density with dimensionless value equal to r 2 . The 
charge separation is of order of the entire film thickness, d. As a consequence, the electric field becomes increasingly 
screened, an effect most pronounced deep in the interior of the system and which spreads toward the two electrodes 
with increasing charge density. This is not difficult to understand. If one groups the layers of opposite charge in 
nesting pairs, with the outermost pair being at z/d = and 1, then the system resembles a set of nested capacitors, 
of alternating polarity, in series. As each capacitor contributes almost no field outside of its plates, the electric field 
can only be large near the two bounding surfaces of the entire system. In the limit of very large densities, it is clear 
that the charge is concentrated in the A lamella nearest the plates, and the field and associated potential drop is only 
significant in the B layers adjacent to the plates. 



B. The Perpendicular Orientation 



For the case in which the lamellae are perpendicular to the plates, we must solve the Poisson-Boltzmann equations 
in two dimensions. Because of periodicity and symmetry, it is sufficient to solve them only in an area of length d in 
the z direction and of length A/2 in the x direction, from the middle of one A lamella to the middle of the adjoining 
B. Equipotential contours are shown in Fig. [7] for an applied potential v — 10 and with a charge density r 2 = 100 
. The boundary between A and B regions runs vertically through the center of the figure with the A region to the 
left, the B to the right. One sees immediately that there is little variation of the electric field from A to B regions. 
This is due to the fact that the lamellae are quite narrow compared to their height. But one also notices that the 
equipotentials are very different from the equally spaced ones which would characterize the potential in the absence of 
ions. Because the top, positive potential, plate attracts the negative mobile counterions, the electric field is screened 
effectively near it and the potential changes slowly; most of the change in potential occurs near the bottom plate. 

As the density of free ions is increased, the A lamellae act more and more like conductors, screening their interiors 
from the external electric field, just as in the parallel orientation. But in the perpendicular orientation, they provide a 
direct path between capacitor plates with the consequences that the external voltage source must place an increasing 
amount of charge on the plates in order to establish the requisite potential difference, and must expend ever more 
energy in order to do so. 



C. The Transition between Orientations 



With the solution of the Poisson-Boltzmann equations in hand for the two configurations, the corresponding free 
energies can be evaluated for any given external voltage, and in particular for voltages at which the perpendicular 
orientation of lamellae becomes of lower free energy than the parallel configuration. We calculate the free energies 
per unit area and /'[ for these two configurations and define 

A/ = / e i - fl (21) 
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Further we note that in the perpendicular configuration Sa — Sb = <5 P iatc because we have assumed that the A and 
B lamellae are of equal thickness. We assume that the surface fields cause the plates to prefer the less-polar polymer 
comprising the B lamellae, as in the experiments of Xu et alP^, so that in the parallel configuration Sa = and 



Sb = 26>pi a t c . From a comparison of the total free energies of Eq. (12) for the two orientations, it follows that the 
perpendicular configuration is favored when 

/ B 2 Xe 2 \ 

A/<- ( 7 x-7fl). (22) 



2c k a 

In Fig. [8] a plot of A/ vs. the dimensionless applied voltage v is presented for the case in which the A lamellae are 
locally neutral. The lamellar periodicity A and film thickness d are taken to be A = d/10, which corresponds to the 
system of Xu et aH Results for four different values of r 2 , the dimensionless density of the mobile ions, are shown. 

The critical value of the external voltage needed to align the lamellae is simply determined by the value of v at 
which a horizontal line drawn at the value A/ = — (f] 2 Xe 2 /2cqK,a)(ia — 1b) intersects the curve Af(v). In particular, 
consider the experiments of Xu et alP3 for which A is PMMA, B is PS, and ja — 1b ~ 0.5 x 10~ 3 Jm~ 2 . With an 
absolute temperature T = 430K, a periodicity A ~ 30 nm and k a = 6.3, the characteristic surface energy 2€qka/ (3 2 Xe 2 
is 5.1 x 10~ 6 Jm -2 leading to a dimensionless surface energy difference of A/ = —97. A horizontal line at that value 
is shown in Fig. [8] We see that the dimensionless voltage needed to bring about alignment is reduced from the value 
without ions of about v = 200, corresponding to Vq ~ 7.5 V, to about v — 60, or Vq ~ 2.2 V, when the density of ions 
corresponds to r 2 = 100. (For a film thickness d = 300 nm, the unit of number density, ka^o/ f3e 2 d 2 = 1.4 x 10 20 m~ 3 
so r 2 = 100 corresponds to p_|_ = 1.4 x 10 22 m~ 3 ). 

This reduction in critical voltage comes from the fact that the counterions do little to lower the electrostatic energy of 
the system when the lamellae are parallel to the substrate as well as locally neutral. In this case, there is no separation 
of charge on the order of the film thickness. But the ions do a great deal to lower the free energy when the lamellae 
are perpendicular to the plates so that there is a separation of charge on the order of the film thickness. Of course, 
this separation of charge decreases the entropy and is therefore opposed by thermal effects. In order to determine 
whether a reorientation will result from the two competing effects, one must carry out the calculation explicitly, as we 
have done, taking into account both the decrease of the free energy due to the reduction in electrostatic energy and 
the increase in free energy due to the decrease in counterion entropy. As the charge density increases without limit, 
the free energy of the perpendicular orientation at any non-zero voltage decreases without limit so that the curves of 
free energy difference approach the ordinate of Fig. 8, and the voltage needed to bring about reorientation approaches 
zero. 

We next consider the less restrictive case when the system is globally neutral. The free energy difference between 
perpendicular and parallel orientations is shown in Fig. [9j Again there is a reduction in the voltage needed to reorient 
the lamellae, from v ~ 200 (r 2 — 0, or no ions) toux 160 for r 2 = 100. This is a much smaller reduction than in the 
locally neutral case for which the reorienting potential was driven down to v 60 for r 2 = 100. The reason that there 
is any reduction in the reorienting voltage in the globally neutral case is subtle for there is a macroscopic separation of 
charge in both orientations. However, while the charge in the perpendicular alignment can be separated by a distance 
d, in the parallel alignment it can only be separated by a smaller distance d — A/2 when the B layer is next to the 
plates. Hence, the ions lower the free energy of the perpendicular configuration more than that of the parallel one 
with a consequent reduction in critical voltage. Because the difference between charge separations of d and d — A/2 
is relatively larger the smaller the value of d, we expect that the relative reduction of the potential needed to reorient 
the film is greater for thinner films. We have verified this conclusion by calculating results for films in which d was 
reduced to 4A. Whereas in the d — 10A film there is a reduction in the reorientation potential of about 17% when the 
dimensionless surface field is equal to 100 and the density of ions is r 2 = 100, in the thinner film the reduction for the 
same surface field and ion density is about 32%. 

Were the A material preferentially adsorbed to the plates, there would be a macroscopic separation of charge of 
order d in both orientations. In fact, the polarization would be somewhat larger in the parallel orientation in which 
the free charges could coat the plates entirely as contrasted with the perpendicular configuration in which they could 
coat only half the area of the plates. Thus we would expect that the addition of ions actually enhances the parallel 
orientation with respect to the perpendicular and would cause the critical voltage to increase rather than decrease. 
We have verified that this is indeed the case, although the increase in critical voltage is not large. We also expect the 
effective surface tension of the A block with the substrate to be reduced by an amount cqKaVq / Xd ■ This expression is 
identical to the reduction of the liquid-solid interfacial tension of polar liquids in external potentials (electrowettingj^l. 



IV. CONCLUSIONS 



We have considered a system composed of lamellae with two different dielectric constants and in which one kind, 
the A lamellae, contain mobile negative ions in a smeared continuum background of positive ones. We calculated the 
effect of these ions on the free energy when the lamellae are oriented in the direction of the applied field and compared 
it to that when the lamellae arc oriented parallel to the two external electrodes, which is the orientation favored by 
surface interactions. We have assumed that the system is in thermal equilibrium so that the comparison of these free 
energies determines the magnitude of the applied voltage necessary to bring about the reorientation of the lamellae. 

The effect of the ions can be quite significant, but we have shown that it depends on several factors. The ions have 
their largest effect if the system is locally neutral, i.e. they neutralize the background charge of each of the A lamellae 
of the film. In that case, a macroscopic separation of charge, with its concomitant lowering of energy, can only be 
brought about if the lamellae undergo the desired reorientation. We find that for the surface interactions commonly 
encountered, the voltage needed to bring about the reorientation could be reduced by more than half. These results 
are in fine qualitative agreement with the experiments of Kohn et alf^ whose system was one of local neutrality. By 
applying alternating electric fields at various frequencies, they could explicity determine the effect of the free ions as 
contrasted with those originating from the different dielectric constants. They indeed found that at lower fr equen cies 
at which the effects of the free ions become important, the orienting effect of the electric field was enhance d 23 * 27 !. 

For systems in which the mobile ions can move from one A lamella to another so that it is only restricted to be 
globally neutral, the effect of the ions is predicted to be smaller than in the locally neutral case, and to depend on 
other factors as well. In particular, if the system is globally neutral, we predict that the effect is greater if the surface 
interactions prefer the B (no ions) lamellae, for then the separation of charge in the perpendicular orientation is larger 
than that in the parallel one. The relative reduction in the reorienting potential should be the greater the thinner the 
film. If the surface interactions prefer the A layers, then the presence of ions can actually increase the voltage needed 
to bring about the reorientation. 

Insight into the question of whether the lamellae of a particular system are locally or globally neutral can be obtained 
by altering the plates so that they prefer A in one experiment and B in another. A locally neutral system will be 
relatively unaffected by such a change, while a globally neutral one would be significantly changed. An alternative is 
to solubilize a known amount of ions in the pure polymers. Integration of the current that passes through the polymer 
at a given external potential allows one to infer the fraction of immobile ions and the solubilization energy. 

We have noted that separation of charge which occurs on the scale of the film thickness decreases the entropy 
making reorientation more difficult. Hence a reduction of temperature certainly makes it easier to bring about the 



reorientation. In our calculation, this is reflected in the definition of the dimensionless measure r, Eq. (14 1, which 
shows that a decrease in temperature is equivalent to an increase in ionic charge density. 

At constant temperature, all properties of the system depend upon the ratio of the ion density to the applied 



potential, p + /Vq, as seen, from Eq. (17 1. Hence for a system with a given ion density at a fixed temperature, any 
"effective" dielectric constant must depend upon the applied voltage itself, and is therefore not a particularly useful 
concept. 

Finally the effects of the ions in our system cannot be attributed to a change in the strength of the interaction 
between components, so whatever role such an effect may play in diblock copolymer systems 1 ^, it is supplementary 
to the equilibrium effects considered here. 

In sum, we have calculated the change in the aligning potential that can be expected upon the introduction of ions 
into a patterned system of different dielectric material^ such as block copolymer, and have highlighted some of the 
issues to be addressed in order to exploit such systems for further applications. 
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Science Foundation under Grants No. 284/05 (YT) and No. 231/08 (DA), and by the German Israeli Foundation 
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Appendix I: Derivation of the Electrostatic Energy 



To obtain the partition function of the system, one must first obtain the energy of a given configuration of charges. 
The electric energy is specified by the number density of the A_ mobile negative charges 



JV_ 



p_(r) =^<5(r-rj 



(LI) 



where is the position vector of the i'th charge. The energy is obtained in the standard way and leads to the usual 
result for a linear dielectric 



U 



\e Q J K (r)[W(r)] 2 dr. 



where the potential, V[p_], is a functional of the number density via the Poisson equation 

1 



V7(r) 



47TK(r)e 



p_(r). 



(1.2) 



(1.3) 



This energy, however, does not include the work that the external electric power source, which is part of the system, 
must do in order to keep constant the potential difference between the two electrodes^. We include this contribution 
to the energy and denote the result for the total energy of the system Utot- Using the fact that one plate is at potential 
zero, and the other at potential Vo, the decrease in total energy can be written 



Utot-U = - J V a(T\\)dr\ { , 

= J V(r,|).D,(r||)dr|| 1 



(1.4) 
(1.5) 
(1.6) 
(1.7) 
(1.8) 

In Eq. (II. 41) above, the charge density of the upper plate is denoted cr(rii ) with rii the position vector in the plane 
of the plate. The Maxwell equation V • D = epf ree has been used to express this charge density in terms of the z 
component of the displacement field just outside the plate, D z , which appears in Eq. (1.5 1, and again to relate the 
divergence of the displacement field in the bulk dielectric, Eq. (1.7), to the free charge there, Eq. (1.8). The second 
integral in that equation is over the regions of A lamellae which contain the charges. Combining the expressions of 
Eqs. (1.2) and (1.8) we obtain for the energy of the system the expression 



J V • [V(r)D(r)]dr, 

J [W(r)] • D(r)dr + J V(r)V • D(r)dr, 

e K(r)[W(r)] 2 dr+ / V(r)e[p + - p_(r)]dr 



Utot[p-] = -\ J e K(r)[W(r)] 2 dr + J V(r)e[p+ - Mr)]dr 



(1.9) 



which, like the potential V[/3_], is a functional of the charge density p_(r). The partition function in the canonical 
ensemble is 



Z(T,n,V ,N + ,S p i ate ) 



AUA 



Kn7 I njlidr^ exp[-l3Utot[p-}} 

Ad J 



(1.10) 



where As is the de Broglie wavelength of the negative charges. 

There are many methods to obtain from the exact partition function the mean-field approximation to the free 
energy, such as that pioneered by Edwards^ and used commonly since^HH They all lead to the result 



F el (T,fi,V ,iV + ,<W) - U tot \p.(r)]-TS, 

= -\e ( n{v)[VV{v)] 2 dv + e I V(v)[ P+ - p_{v)]dv 



k B T [ p_(r)]n[p.(r)/p + ]dT, 

J A 



(1.11) 



(1.12) 
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where p~(r), in contrast to p-(r) of Eq. (I.ll, is an ensemble-average density 

P- (r) = P+» f , L ' 'I r m A 

J A drexp[j3eV(r)\ 

and is the source of the potential via the Poisson equation 



(1.13) 



V 2 l/(r) 



47TK(r)e 



(1.14) 



The mean-field free energy, Eq. (1.12), is seen to be simply the total electrostatic energy, Eq. (1.9), but as a function 
of this average charge density, supplemented by the contribution to the entropy of the mobile counterions. 
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Appendix II: Numerical Solutions 



The rescaled electrostatic free energy / G i defined in Eq. ( 11 ) can be written in terms of the rescaled potential W by 
using the relation p-(r) = p + 9 cxp(W(r)) to eliminate the explicit dependence on p_(r) in Eq. ([£]). The result, after 
changing to rescaled variables and neglecting a constant which contributes equally to the parallel and perpendicular 
configurations, is 



1 r A ' za r L k(t ?\ 

/el,global = - x / / -^^{VW(x,S)} 2 dxd~Z 



r W(x,z)dxdz 



A 

+ H^) ln * <ni) 

This expression is valid for both the globally neutral parallel configuration and the perpendicular configuration. In 
the case of locally neutral parallel lamellae, each lamella labeled by i has its own constant 9i and the rescaled free 
energy is 

\/2d 



i r A i za r L k(t ?) 

/eUocal = - q / / ^^[VW{xrz)] 2 dxd~Z 

£ Jx=0 Jz=0 K A 

+ r 2 J J W(x,z)dxdz 

where Qa,% is the volume taken up by the ith lamella within the integration region f2. These rescaled free energies can 
be calculated once the rescaled potential W and the constant 9 (or constants Oi in the case of locally neutral lamellae) 
have been determined by numerically solving discretized versions of the equations 

\/ 2 W{x, z) = - r 2 (l - 9 t e w{i ^ in lamella A* 

V 2 W(x, z) = in B regions (II.3) 

subject to appropriate boundary conditions. In the case of parallel lamellae where W is a function of z only, the 
following boundary conditions hold in both the locally and globally neutral cases: 

W(0) = (II.4) 

W(l) = v (II.5) 

KTW'fr)^ = K + W>(z t )l ight (11.6) 

W(zi)U = W(zi)\ lisU , (11.7) 

where the AB interfaces are indexed by i and and k + are the relative dielectric constants on the two sides of the 
interface at Zi. In the case of globally neutral lamellae, the unknown 9 is determined by the boundary condition 

W'(0) = W'(l) (11.8) 

In the case of locally neutral lamellae, each unknown 9i is determined by a boundary condition describing the neutrality 
of the zth A-lamella: 

W'{cn) = W'Qh), l<i<n, (II.9) 
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where n is the number of A-lamellae and dj and hi are the scaled positions of the edges of the ith. such lamella. For 
the perpendicular configuration, the boundary conditions are 



K A 



dW 

dx 



dW 



dx 

dW 
dx 
dW 

dx 



=A/4d+ 



x=0 



0, 



W(X/4d,z)\ ieft 



=X/2d 

W(X/4d,z)\ lisht , 



-k b 



KA L { 


dW 

dz 


2=1 


i-A/ad / Qw 




dW 


Ix/id V d ~ z 


z=l 


dz 



W(x,0) = 
W(x,l) =v 

< z < 1 
< z < 1 



< z < 1 
< z < 1 



z=0 



2=0 
f/.f 



0, 



(11.10) 

(11.11) 

(11.12) 
(11.13) 

(11.14) 
(11.15) 

(11.16) 



where x = is in the middle of an A lamella, so that x = \/Ad is the location of the AB interface. Eqs. (11.13) 
and ( II.14[ ) enforce the symmetry of the system with respect to reflection about the center of each A or B lamella, 
respectively. Eq. (11.16) determines the value of 9 by enforcing the neutrality of the system. 
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Figure Captions 

Figure 1. Parallel case, no ions (r 2 = 0). The dimcnsionless applied potential is v = 10, and k = ka/kb = 2. 
The sample thickness and wavelength are related by d = 4A. The dimcnsionless potential, W = jdeV is plotted 
vs. the dimcnsionless coordinate z/d. A- type regions are shown in white, B regions in gray. 

Figure 2. Parallel case, local neutrality. The dimensionless applied potential is v = 10, and k = ka/kb = 2. 
The sample thickness and wavelength are related by d = 4A. The dimensionless electric field, —W, is plotted 
vs. the dimensionless coordinate z/d. A-typc regions are shown in white, B regions in gray. The case for no 
ions, r 2 = is shown with solid lines, while r 2 = 100 is shown with a dotted line. 

Figure 3. Parallel case, local neutrality. The dimensionless total charge density, —W", is plotted for v = 10, 
k = 2. The sample thickness and wavelength are related by d — 4A. Solid line, r 2 = 50; dashed line, r 2 = 100. 

Figure 4. Parallel case, global neutrality. The dimensionless potential W for v — 10, k = 2, and with r 2 
varied. Solid line, r 2 — 0.01; dashed line, r 2 — 50; dotted line, r 2 = 100. The sample thickness and wavelength 
are related by d = 4A. 

Figure 5. Parallel case, global neutrality. The dimensionless electric field, —W, is plotted for v = 10, k = 2, 
and with r 2 varied. Solid line, r 2 = 0.01; dashed line, r 2 — 100. The sample thickness and wavelength are 
related by d = 4A. 

Figure 6. Parallel case, global neutrality. The dimcnsionless charge density, —W", is plotted for v = 10, 
k = 2, and with r 2 varied. Solid line, r 2 = 50; dotted line, r 2 = 100. The sample thickness and wavelength arc 
related by d = 4A. 

Figure 7. Perpendicular case. Contours of constant dimensionless potential, W, at intervals of 1.0 are plotted 
in the x/d,z/d plane for v = 10 and r 2 = 100. As A = d/4, we need only plot x/d from 0.0 to A/2c? = 0.125. 
The interface between A and B lamellae is at x/d = 0.0625 with the A region to the left and B to the right of 
it. Areas between equipotentials are filled with different shades of gray for clarity. 

Figure 8. The difference between the dimensionless electrostatic free energies per unit area (without surface 
fields) of the perpendicular orientation and the locally neutral parallel orientation is shown as a function of 
dimensionless applied voltage v for four values of ion number density; r 2 = 0.01 (solid line), 50 (dashed line), 
100 (dotted line), and 1000 (dotted-dashed line). The sample thickness and wavelength are related by d = 10A. 
The horizontal line is at the value A/ = —97 (see text). 

Figure 9. The difference between the dimcnsionless electrostatic free energies per unit area (without surface 
fields) of the perpendicular orientation and the globally neutral parallel orientation is shown as a function of 
dimensionless applied voltage v for four values of ion number density; r 2 = 0.01, (solid line) 50 (dashed line), 
100 (dotted line), and 1000 (dotted-dashed line). The sample thickness and wavelength are related by d = 10A. 
The horizontal line is at the value A/ = —97 (see text). 
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